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ABS TRACT 


Various aspects of reliability test programs when the number of 
1tetmS tovbe tested equals or exceeds the mumber of test chamber 
available for testing and wmen the underlying failure distribution is 
exponential are studied. Specifically, two bivariate processes are 
formed and characterized; (i) the time until the r“ failure and the 
accumulated test time of the first r failures and (i1) the number of 
failures in the time interval (O,T) and the total ese time accumulated 
in(O,T). These are treated in the two cases; replacement andii@an 
replacement. Relations which exist between contractual specifi- 
cations for reliability demonstration and the underlying distributions 
of the test program are discussed. Maximum likelihood estimators 


are developed for the distributions when necessary . 
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I. INTRODUCTION 


A reliability program manager shall be considered as one who 
has the task to demonstrate , according to contractual specifications, 
the reliability of a piece of equipment. Such a manager is commonly 
faced with certain physical constraints such as limited test facilities, 
number of items available for testing, time available to complete the 
test program, etc. The manager is also faced with various costing 
considerations such as the cost of the items to be tested, or keeping 
the limited test facilities in use for a specific equipment vice another, 
etc. In view of the preceeding, the manager is concerned with such 
problems as, "what is the probability of successful demonstration of 
the specified reliability throughout the test program," "what is an 
estimation of the time required to complete the test program," and 
"what is the number of items which will fail during the testing." It 
is the purpose of this paper to provide mathematical answers to 
questions of this type and hence be an aid to such a manager. 

Throughout this paper it will be assumed that all items to be 
tested have an exponential failure density (i.e. f(x) = he hx and \ 
is the failure rate of the item to be tested.) The assumption of 
exponential failure density is often made for electronic devices and 
it has been shown [Ref. 4] that complex equipments tend to have an 
exponential failure density. This assumption also leads to a large 


number of mathematical simplifications. The extension to other 


failure laws is difficult and the present analysis can serve as a 
skeleton for that extension. The mathematical developments in this 
paper heavily rely upon the properties of the Lapeee wransfornmr 

It is a common practice to write contractual specifications for 
reliability demonstration which call for the demonstration of a 
specified lower confidence limit on reliability, R1, with a confidence 
coefficient % where the item is required to operate over a mission 
time T. Lloyd and Lipow [Ref. 4] show that this specification gives 
rise to the demonstration of an upper confidence limit Le for »\ with 
confidence coefficient Y& when the failure density is exponential. 

A 
This means that the estimator n\ for r » the failure rate , must be 
less than Nein order to assert that the specified reliability has been 
demonstrated. 

Truncated life testing [Ref. 4] is a technique of experimental 
design that economizes in terms of either the number of items 
tested or the amount of time spent testing, or both. This paper is 
concerned with those designs that terminate er either a 

(1) fixed number of failures or 

(2) fixed amount of time. 

Other (sequential) stopping rules are popular [Ref. 4], but will not be 
considered. The design of the experiment will trade cost against the 
utility of Success. The total cost is the sum of the number of items 
tested multiplied by the cost per item and the amount of time spent 


testing multiplied by the cost per unit time (for use of the test 





facilities.) Thus the total cost is a random number because the time 
spent testing for stopping rule (1) is a random variable and the 
number of items tested is random for stopping rule (2). Thus the 
total expected cost is a function of » and the number of failures for 
stopping rule (1), and is a function of N and the termination time for 
stopping rule (2). Generally the expected cost will increase as the 
failure rate mn decreases. On the other hand, the lower the failure 
rate the more desirable is the equipment and the more likely it is 
that the contractual specification c\e do can be demonstrated. It 
follows that the joint distribution of (1) the estimator and the terminal 
time and (2) the estimator and the number of failures will play a key 
role in balancing cost of testing against the value of successful 
completion. This latter value is a utility type consideration and 
must be developed separately for each particular kind of item. A 
Similar statement applies to the cost coefficients. It should be noted 
that many modern systems have high item cost and high test facility 
usage costs. Only the methodology for treating the joint distributions 
mentioned above can be treated in general. Also, the program 
manager has the option of terminating a test program early if the 
chances of Successful completion diminish during the course of 
testing. Such a decision should be based on probability statements 
concerning the aforementioned random vectors conditioned by the 
current status of affairs. Thus an understanding of these quantities 


is needed when they are viewed as bivariate stochastic processes. 


It is assumed that the program manager is constrained by a 
limited number, c, of test chambers. If the number of items 
available for testing, i, is less thanc, then all items are put on 
test at the same time and testing takes place without replacement. 
Without loss of generality this case can be indexed by i=c. 
Otherwise i>c and the procedure is to test with replacement until 
the reservoir of Spares is depleted and then continue to test without 
replacement. The case i=cis treated in Chapter II and the case 
i>c is treated in Chapter III. The joint distributions and the 
marginal distributions and the moments for each type stopping rule 
are developed. It has been shown by Read [Ref. 5] that the normal 
approximations to these distributions are valid if the number of items 
failed and the failure rate are not too small. Explicit distributions 
will be necessary when these quantities are small. An indication 


for the necessary size of these quantities for the normal 


approximations to be valid is also given by Read [Ref. 5]. 
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WW. “hie WiMee OF TTEMS TO BE TESTED 


DOE S' N@T EX@EBIEID THE NUMBER GF TES W CHAMBERS 


me. SMOPPING RULE! STOP Fadk URES 

In this first case (i=c or p=0), all items are put on test at the 
start of the testing program and the test is complete when r of the 
c items have failed. 

Liot 4 = Tiere to'first fabilume from t=O 


fog Time to second failure from t=O 


Tp, = Time to rth failure from t=0 


Srp = Tp 7 Tr 


then S. has density (c-j+1) , e 


-(c-jt1) dx, 
j nlm shee > 


x > 0. The Sjrane 
independent and hence the Laplace Transform of the joint density of 


the S;'s is: 


r 
c-jt1) A 


Ptisyessy- - -, SA) = 3 3 
j=1 sj + (c-jt1) 


(1) 
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It is of general interest to characterize the distribution of the 


order Statistic. Jbepipegieaolace | ransfOpmuon the joint density of 


T4> Toserss ile. be: 
» 
ae uy Ty 
£ * Cees = £ Le pons 
J 
sco Ty 2 ey , and hence, 
a= 


ie 
_[ -2sle 6) 
f* (ua va = oa e. j=! a: 


= £ [Ele | 


% 
Therefore , letting 4) = 2 Um yields: 
=) 


£¥ us uyedy) = Tr Co~jts) A 


- (2) 
“i (c= ja\d + 24% 


ky 


Inverting equation (2) on U4 gives: 


e 
= { (cry) +Ur + 2. mel T 


Ge We Ee CoD ee 


eat { (cA + Ur 4 J are} 


K=% 


ns 
(c~j+t) A 
ip a 


J23 (eeges\At 2 Uk 
k=) 


NZ 





and noting that the Laplace Transform of 


~ AX 
€ U (xe Eis) 


where (U"(:xrty ) is the unit step function translated to - ty ,is 


es (A+) CT, 


= 





(aA+~\ 


gives rise to the joint density of T,; To> eee, T, as: 


ay T eyeremvans iy 


ee or vit Uy = i : : 
Ae \ eae ee te (3) 
aon 





for O ¢ty (to Ceoe Cth 


th failure can 


The marginal distribution of the time to the j 
be calculated easily from equation (2). 


It has been shown [Ref.2] that the maximum likelihood estimator 


for A upon completion of the test program is: 


A 
Nh = (4) 
where B. is the total accumulated test time of all items. 
(Alternatively, this may be shown from equation (3)). Thus the joint 
NN 
distribution of A and the time to completion Ty, is tantamount to the 
joint distribution of T, and B.. This joint distribution has been 


characterized by Read [ Ref.5] and for completeness the Laplace 


Transform and the density is recorded: 
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~ Ce yethd 
f(u,u)= || we Oma (5) 
j=) 


u r-\-) -X, A Y-2 
Fr bux) s oy 2. (1) : (;') Ce [ Ae (c-3)] : (GF 
. (y-2)! 





Uy. - x, (e-g)] 


where U(x) is the unit step function. 

When the program manager desires to mark progress at some 
time Taurine the course of testing, the data available at that time 
will be a complete history of the testing including N-yv and By , i.e. 
the current tally on the number of failures and the accumulated total 
test time. Due to the fact that the exponential distribution has no 
memory, all data other than N +y and Bry are irrelevant in terms 
of making updated probability statements concerning T,, and B, - 

It also follows that the conditional density of Ts, Br given that 

Pine =m (for m<r) and By =b can be characterized from equations 
(5) and (6) by replacing r with (r-m), ¢c with (c~m) and by translating 
ec ea late etsemiarrees ir iols  iilalehee 


. * (7) 
Teac Ceca Cee, Hae —_ 


a Teurkh wm 


e (c-m-j3 +I A 


ey Mt (Unt n) Co-m-j et) 


and letting K =r-m yields: 
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tr By | Nw =m, Beeb Oxi, ¥.) > 


is~) | 7 (X,~lb) 
\< paneer Keo) : 
Kd” ( oye 4 («-') e 8) 
K+ 2 
TL On-b) - (u-T)c- m-y)1 U fa = 8) - (k,- 7) (omy) 
Ck> oe 


The univariate distributions developed in the Same manner 


from equation (5) are: 


v7] 


=i Coe ~ AX 
fru) =[s)r re Lire ] je @& 
The |Nqem, Beeb Oh) = 
one (xX, TV(e-r +1) SX os KK-] 
ape mn © Lie ] m7 ee 
\< 
vol aXe 
_ AOR) 
sree Core = — ae \ Xue7ee@ (11) 
Whee ee oye 
dhe ay 
alll ’ Xai (12) 


(\<~1)! 


The moments needed for mean value analysis and for the normal 


approximations developed by Read [Ref.5] are: 


NV 
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2 = 
ET] Ny=m, By=T} aes A 2- (Cs jenah) 


‘a 


= | / 
VaR i Ty] Nezm. Be = | = — Se 
{ r v ‘aes } \> jem) (o-y +) 


£{ Br} = — 


VaR { Br} = re 
/ 


Val Ca 


d 





E{Be | nypem, By= THe TH 





Var { Bel Neem BeeT} = YM 
\> 
y=l 


Cov) Tr, Be} = A 2. __ 


Ne ayo Te! Gi) 


yom! | 


Cov {Th Br | Ny=m, Byes Tt = c. ee 


\e ae Ce-mny ) 


(13) 


Now consider those problems posed for the reliability program 
manager. The manager can successfully demonstrate reliability if 
A A 
he can produce an estimator for \ , say \ , such that A £ Aco 
where Ag is an upper confidence limit for x. » with confidence 


coefficient ¥ . Now upon completion of the test program, the 


manager will have successfully demonstrated reliability if: 
. r 
=a (14) 


Therefore, it is required that the total test time accumulated during 
ane test program B,. must be greater thanir7 a, were B mastdcissne, 
given by equation (11). Hence, if the program manager is willing to 
use a design estimate, say op Pmole » he has available to him, 
the probability of successfully demonstrating reliability, an estimate 
of the total time for the test program, and various other pieces of 
information prior to the start of the testing. 

Now given that at some time, 7 , prior to completion of the test, 
the manager desires an estimator for \ based upon the data which is 
available at ‘y , it has been shown [Ref. 2], that the maximum 
likelihood estimator ie is 

cae 


TO Ba 


A 
and by simply upgrading the distributions in question with \y P 
the program manager has a continuing estimate of the most important 


characteristics of his test program. 
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B. STOPPING RULE: STOP AT A FIXED TIME T 
The distribution of the order statistic T,, To, oe Le given 


N-- =r has general interest and can be obtained from equation (3) as 


follows: 


fay, ere TANS Oe Oe oN oe 


Prob | N+ z+ lens > Tr} TT (Tr Trey Te) 
Prob {N+ =v} 


Noting: 


aac.) Ae Ty) 


‘Prob { Nrer Pare = ty} = = ’ (Ge as 
a 2) eee 
Prob { Nrer } = es (| - @ es e 
and therefore 
ile | Nyr=r C0), be, ey Te) S 
? (15) 
Se 
a) ieee 
i-. 
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Under a fixed time stopping rule, with 1=c, it has been shown 
(Ref.2] that the maximum likelihood estimator for \ upon completion 
of the test program is: 

Ar 


‘\ , 
(16) 
B+ 


and the joint distribution of N-; B- is needed for experimental 
G@esign and decision making. Recognizing that Nai ws a discrete 
random variable while B— is a continuous random variable, their 


+ 


foie Censpby will be referred to as; 


pi (Prob (Neer, By éb)) 
F ? ( (6 eee (17) 
TRF aa 3 b 


To develop this quantity, the following relationship is needed: 


B- = By, +(T-T C-N 
i; Eo Nee Tt? 


and hence: 
= (N+ =- : B+ & b = 


ia ( lis foul ry ee Sie ee (c-r) = b } 
letting J. = B. + (CT — T,) (c-r), equation (17) can be written as 


Prob tr £T, T > ee. é b) 


mal 
and this probability is equal to: 


b 


if fore) 
\ dn, | cli \ Cl Xie Kail ) ale 


© 


is, 


where f(X45 Xo, X) is the joint density of T,, A ase” ene 
let the Laplace Transform of this density be: 


eens muy beat — es [ Jv ve eon 


: ) (18) 


Se Se 
Hea, u,\ = € i 


and noting the following relations: 


r 
VE = ‘y S, 
a 
v +) 
Tra) = S, 
w=! 


de - (6-r\ 7 = va (r=) 41) S\ 


! 


Jz 


and hence make the following transformation of variables in equation 


(1): 
ees ee 


Ay = Aye Uy + Uy (r-y 4) 


Pree F Ue 
gives the Laplace Transform of f(x4, Xo; Xn ) as: 


~ Uy (C-r) im 


Cur wy, ay) = a (c-r)\ “St y 
(c-w) (ug + Ae a. 


v-{| (Bi) 


a 
Us + qc-r) A ps, Ceo aie (u, + Upt Uy (v-y) + CO-) NA) 
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Now inverting on U4 and denoting the result ly, f*, yields: 


m - Wa (c-r) T 
Ls, {* = \* (c-r) ie) € . 
(us + a (uz 4(ce-r)A) 


-1 ae - 
2 _ ( 4 ) ’ 


= Ut 
Now by factoring e igh ba Er Neiee from under the 





~ Xx, Cur t*aagl(vryy + (c-ysd J 


summation sign, and applying the binomial theorem to what is left, 


yields: 
“ “Ux le-r)T  _¥j gr _—K Ue 
lu,f* = A "T (e-r) a%) ©. ee 2 
(U3 + yy 5° [ uy + ce-v) 4 


[ ot htarh) - | ¥=| 


Next invert on Uo: 





lu,u, f* = 
—uU Ce-r) 1 -K% Gav - xX Cc-v) J 
Na e-ey < (<) > Sia i) 
Y 
Tro Va 
~Xr X, (Uy +d) ie 
e a - : ~| } U (Xi - X,) 


~ Li, (eer) T ~ Xv (err) ) 
e = 


¢ 


w+ | , 
- & tee) £< ) Seles 
(us +A) 


~~ X\ + = Cu «d\ r- | 
e ’ “ " ute eS . a U (X2-X,) 


a1 


Hence 


T fore) 
\ \ Tous ae ip dX. dx, = 


0 a 


~ ie (dana? oa TRUER ei 
e Us (C 5 Cc" ( 3t } (19) 


r 





(s 
( lA i yO 
which can be rewritten as: 


n’(8) 5 (-))9 («) e- T Ce-j)Cuy+d) (20) 
jr0 





(Uy +n)" 


Therefore , upon performing the final inversion: 


Ly. cael ei (S) 


Yr 


ve i ' C~ \ 
> (-iy 4 { *.) Lb-Te-yy) U (e-TKe-y)) <1 
Ve Lv <1)! 
bead 
which 1S the desired joint density. 
Now by considering the fact that 


re’ Neer} = (Sa aes = Ast ae 


and dividing this probability into equation (19), yields: 


-T (ug +r) 


une = _ Us (C= 
Eje ° "| Nrer| = e | MWe 3 bc - 


iy + d 
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and noting 


-2. Np ~ux,B 
E{e* = 


Me 


a ee { Nr] g{eo" Nr] 
@) 


C SC ee ( 
7 : (¢) \ | eel Nor _ Uy fe" Ustr) “hil (22) 
= Cus tro 





which can be rewritten as: 


~ T (ugar “wr TT yc 
e prove heme a) -1] (23) 
ere ete SP ) 
(us 4A) 
The first and second moments are obtained from equation (22) 


aS: 


Efvrp = ¢ fi- en" J 


e{B}=s ct - & fat leer 4 


(24) 
Van {Net = ce ie eT) 

Van {Br} - >, 4 oa. ean 

Gans [Nr Brl ' a ~~ . OT a) 
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For purposes of marking progress when the elapsed time of the 
tesulime isere |, tie memeevyless property oF Enevexpomential 
distribution again makes the updating of probability statements 


relatively easy. Given N =m, B ~~ = Do, the distribution of 


- 


N+ and B_ can be Obtaimed by replacing r by r-m, c by cm, Ti wits 


i, 
T —% and translating N-_y~ =ro-m, Boley =x — bo- 

At this point the reliability program manager is again armed with 
all the instrumentation required to have continuing estimates of the 
most important characteristics of the test program. It should be 
noted for this case however, that the estimator for A » say \ » is 
now the ratio of two random variables and hence it is this ratio which 
must be less than the A Q designated in the contractual specifications, 


teen? 


N+ 


Ao 


IX 
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ii. THE NUMBER OF eaMwiSene ee tesSTreo 


EXCEEDS THE NUMBER OF TEST CHAMBERS AVAILABLE 


=~. SSOPPING RULE: STOP Ab FAILURES 


The development in this chapter will generally follow the 


development of Chapter II. Testing takes place with replacement as 


long as possible and testing is complete at the rth failure. At the 


start of the testing, there are i items available for testing but only 


c test chambers available, where i > c. This chapter is restricted 


to the case where r > p=i-c. To start, let: 


and 


then Sj 


mas density (ctp-jt1) Ne 


Ty = Time to first failure from t =O 


To = Time to second failure from t = 0 


Tr = Time to al failure from t =O 


So=T,-T, 
Sp =T,, eV n=4 
SMX 
has density c Ae, X>Oforj=i, 2, ..., pand 


—(cto-it1 A x 
Pan ,x > Ofor j =ot1,...,r. 


Hence the Laplace Transform of the joint density of the S,'s is: 
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i cA (c+p-sti)a 


Cr eee) oo a . —_—. ao 
a2 se UT Sit CA s=Pt! Sy epee rt! JA 


Therefore the Laplace Transform of the joint density of 


(25) 


a 


ick, Tose++5 1, becomes: 


P A r (crept) A 
c ana - J 
“ees = || ——= eee | : (26) 
Fl cht Dv, jepu CPG eye ly. 
k= 4 eat 
J k=) 


Now inverting equation (26) on uy yields: 


p (cA + Uy + he ‘i 
al on 
1-3 k= 


r 


(c+p-st+i) A 


; a 
faethe) tidied 1, 
key 


Hence the inversion on U4> Has eee ue produces: 


eX 
L (Esta, + + tp Upgyye Ey) = 


M1) Ep . (c+e-y+i)a 


feper (O+P-S+I)A + y Up 
K=q 


(ay e 


andO < t; < ty oe ee o 


and the format has now become that of equation (2) and 


memee thegietmecensity of 7 


Me e 


T x 5) as eee S 

ey mJ 

-\ 2D tf -(crepa)at, 
€. e A 


mF 
F (tyty- ss tr) = coi 
(cor 7 P)! 


j= etl 27) 
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The maximum likelihood estimator for \ Using a stopping rule 
based on r failures and assuming r > pis developed in Appendix B 


aS: 





> 


Thus again the joint distribution of and the time to completion is 
tantamount to the joint distribution of Tp, B.- 
The development of this chapter differs slightly from that of 
@mapter Il in that the Laplacesiransform of the jotmedensimy of 1,,; 
ae Tant» B,» and J (to be defined later) is developed and from > 

this Laplace Transform the desired combinations of distributions 
are extracted. This method allows for a clear presentation of the 
relations existing between the various random variables involved in 


the process. 


Consider the joint distribution of T,, B., N Bw 


CP? ae 

and again note: 

Prob (1, 4t, B.£b,B.4z,N,=m) = 
Sint | Umer th eee een Se is ule) z | 
where J = Baa CG ame (i-mM) forse < om ete 
(The case m ~ p is straightforward and will not be treated 


explicitely.) 


This probability is equal to: 


Tt v oo b Z 
\ dx, \ dn, | dns | dx, \ i F(X, s%29%3 1 %q aXe) 
oO - fo) O 


@) Ae 


2/ 


wnere F(X4,6+. »X5) is*me joint densityor Lise hes ae 4 i= he 


Now let the Laplace Transform of this density be: 


* — 
Z (Uy see Ue) = 


~ Us (aem) YT 


e @ 


(vl 


MU Te Ur T= Us Tima = Uy Be de (Simm Gm) % (28) 


and note the following relations: 


is 
TO. = =: 
er j=1 J 
TW = = 
= 1 eee 
taat 
ar = Pa > 
E ie 
B. =cCc ? =. 4 > (Camila) =e 
all j=pr1 
—P mM 
J (ism) “CS =(mep * Sy a '” (Cmrjti) «Se 
m ae) a J 
j=l j=pt1 
Hence the transformation of variables 
= + + + + — i= 1 eee 
25 U, Un Un cu, cm PJu, (=. sD 
Sj uy ale Ls a u,, a (e-p-jttu, (tir joa U . jae", .. > 5am 
S}) = uy tug t (ctp—m)u, j=mt+1 
Sj, =U, + (ctpj+1 uy jemte 2... ie 


28 


can be applied to equation (28) and (noting that ctp =i), the Laplace 
Transform of equation (28) becomes: 
f*(uy,+++,U5) = 
i 
Tee t 
= O-( boat cA aa 


e ——— ae ee 
U, + U, + Ug + (U,+A)C + (m-p) Us 


TT (L-s41) A 
J p42 U,+U,4 U, 4 (Uy4d) (b-j41) +U5 (m-pti) 
(t-m)A CB 
TT Gaon apy, 
Jjemrr ect (vy4A) (i-rTtt) 
Now setting U5 - ers = ijee Ue = O yields the Laplace Transform 
of density of dara ‘to be: 
Ptrl ‘ia. 
x ” 
f (44 0/0,0) =(—£5 | vy Attlee (30) 
Ur+ ca Dt tate pe 


Jeptr 
Equation (30) shows that T,, has an Erlang distribution with r 
stages and the indicated rates. Now equation (30) can be written as: 


foe. , O, 0, 0, 0)= 





(31) 
(l-e-1) | a , 
-™ Gia pt ' aS coe = 
fpr (uted) (Ut (tp +1) AJ 
4 
Where Cis 4596. 33 
J Ca ee ee 


29 


and using theorem 1 of Appendix A, the density of Tp, becomes: 





(t) fea i)) = |ftt — a 
es Uke apt Las ( e+ | + 
| je pte 5-p-!) 
(32) 
-cAt p | >-k  fek peak 
e- ies) 7) } 
eae pet) 
(p-5+1) ie (e-k) | 


>imilarly, the settiiag of UR= U5 = Ug = Us = O in equation (29) 


yields the Laplace Wiranste@maieimmeneecemsity Of B. torbe: 





% A ‘a 
f (0,0, 0, Uhr 2) = ( a 
Cieat a (33) 
which is the Laplace Transform of a special Erlang distribution 
with r stages and rate A and hence the density of B, is: 
r-l - AX 
x) e 
fix) = AlaAx) e , x > 
ane (34) 
Now considering equation (29) with U5 = Ug = Ue = O, the 
Paplace Mransiogan Of the joint cdemsity of |, a iS produced: 
F™(U4s O, O, UaA>s 0) - 
! 
sete Hirt aie a 
35 
Ut (Uyta)c ae UW, + Cenpot) ( Uaad ) se 
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which can be rewritten as: 


¥ oc lire)! 
Cae | (an UHR cde (uy + ies 
v-\ 


S en ae | | 


Wier -u mereyen  [ clk, (us Se | 


Sept 
maiererore: 
Iuy f* = (36) 

ie P+ | , ile ao me 
N ¢ (ASPs We | > ole 
ae. (uy ery O7P® Jtpri (j-prat tyy-): 
e ~ X, | Cuyray (a-y)] 
+> 


(uy try"! Cy-pyP™ 


~X Clr Le 
———— 


of S (-1)* } X (uy +d) (P ay 


Cee yh Cae (ae " 
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Now by inverting onuy, — the joint density of T,, B, is obtained as: 
om 


lu, uy, fe Golan) = nee, (37) 
. ‘a — r ' 
eae een.) § > i ] 
i : (L- ry) ei) (jae ; 
5 fept jor eyo! (ra 
( t r-2 : t 
ped ( x -x, (l-f) } LJ (eee xitnaaieds) 





(r-2)lo (f-p)r™ 


i ; ae 
Sent etek Ota PT ur eae 
kKz=o0 (r-p+k-2)! (p-k)] Cp-4) <4! 


By considering the logarithm of the transform, Read [Ref .5] 


has shown: 














9 | roe all 
eee Eh J | . 
A SOS) 
repo t 
Tl ge | (38) 
(cA)? a a A 
. _ ~ 
eee 
Var “ce a 
A rs ! 
| = Loe _ 2 | 
eisupile ly) Sm) oo At Dye ¥) (c-5) 
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The program manager iS now armed with the required 
distribution theory to determine if the reliability testing has 
successfully demonstrated the specified contractual reliability. If 
the manager desires to mark progress prior to the completion of the 
testing, the data-available at sometime , eit » will be a complete 
history of the testing includingN~ , Bey . Again the memory- 
less property of the exponential distribution makes N~q , Bry the 
only relevant data in terms of making probability statements 
concerning the density of T,, B,- Appendix B develops the maximum 
likelihood estimator for » say 7 , for a fixed time T , as: 

7 = iio 
Bage 

It also follows that the conditional density of Tr, B, given 
Nee = m (form & r) and B ~ =b can be characterized from equations 
@e) and (34) by replacing r with (r-—m), it with Cizm) and by 
Geemslating T, to T.-ay+ T 5 and B. to B,-p, + 0, keeping track of 


wnether m4 p or m )> p. Thus: 
$ F = 
Th, BEI Nrem Breh Cu, , ey) = 


Se ai a uy b 
C 7, te, Onotun, O)) 


with the appropriate translations on the indices depending on the 


value of m. 
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The conditional moments are given by: 











\ 
_~ P- ‘I 
E {Te | Ny em c<o\} are | 2, (¢-3) 
Y-p-—l | 
cls 
ey, ~ | Nes m Copy} = os ; eo ea) 





VaR {Te | Neem (ep) } = 





v-p-} 
Var ST. | Ny =m(op)t = 7 a -\)* 
Ef Be | Nem, Bree} = b + cae 
VaR {By | Neem, By ='} z = 


Again the manager is confronted with a eal failure rate 
as derived from the contractual specifications on reliability. The 
system or equipment under test in this case must develop a total busy 
time B, which is greater than r/), . At this point the manager is 
equipped with the necessary managerial aids for decision making 


throughout the test program. 
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Be toe inG RULE: =SsiTer Ab AAA eeD TIME. 1, 
Under a fixed time stopping rule, withi > c, the maximum 


likelihood estimator for A upon completion of the test program is: 


a 
4A 


) 


and again the joint distribution of N-, eis required for experimental 
design and decision making. Using the same convention as in Chapter 


II yields the joint density on ae BL a 


3 Prob (Wr=m , By £b) 
jecl 3b 

in this case however the procedure is not quite so straightforward 

as it was in Chapter II. Equation (29) requires some understanding 


of the indices. Proper manipulation of the indices yields: 


c ( Zgane —U, Tmet — Us Tm 
, Us CT } ag | m+ 7 
ig th3t+Ca at san 
where m < p. 
It should be noted that B+ is a degenerate random variable for this 


CaS. 


-Ds Tm ~U, Tm #1 -~Us Je, 


Ele 


: e+! 
ui, a ca (<-1) 4 


rn rn ger 


hes Flees CA U,+ (c-1)) 


Where m=pti1. 
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~ a Im-4& Tiga ~Us Diese 
ig Ve ail ° } = (40) 


— Us CAem)T 
‘< 


\+ 
__ ae 
Uy~tuy t (m-prdug 4+ cA 


TT Coe) N 
y=pre Uz 4 Us 4 Un-yer) Ge F Canmysr\ a 


(A- Mm) X 


Ux + CA-m\A 


where pt1 4 m < i 


The following relationships should be noted 


~ue Br 
F;e | Neer f = Equation (40) with Us ug= 0 


F { eee} = > Prob {Nrvfe{e "| perf 


y= 0 


| . 5 Cora { Nree] are : ene J rev f 
a 
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The marginal distribution of NEE is obtainable by considering for 


Geen r,r=0, 1, 2, --. 1 , themerobani lity N-- =r. For 


r=0, 1, 2, .-. Pp, the process iS a Wersen peegess with rate ; 
BIeince: 
oe (41) 
Prob {Nr - | = CM ae c 
iG: 


a= Oe erates > 


Of ites p + 1 


and from equation (40), the Laplace Transform of this joint density is: 


£* (ug, U3) a ch go Cea a a 
Up + Uy +c Uy + (e-1)h 


Inversion in U5 yields: 


‘ P = Xu eee) 
a, 68 Be gd we 
i [u;*cr | 


Inversion on ug yields: 


= Pit WIP = heh ge Named 
J earres foe wee i (<-1) \ Xz @€ i a 
e 


Sn 


Now: 


P+i . Pel Pr = Koen aay (C=a).A 
CaS gi cee 
ry 
Ox< ye a 
Finally, 
¥ 
ye a Ga -\T Pf (42) 
PANT pert = C © ae > (AT) ; 
f= Oe) | 


Now for rp pti, the Laplace Transform of the joint density of 


Aner. ale is given by equation (40) as: 


Pr? el 


P+\ 


Uy oe @) 


$* Cur, us — C d 
Av + 


‘a (43) 


TT W-\ +1) 


ie Us + Uy + CanmyrlbA 





Li-~vdA 


Ux + (A-" YX 
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Rewriting gives: 


C 


uF 


Riel 2 ee reer ya | ae | 
colon? 
Ursa ee lber |) 7 ( F-p-2) } (r-7) 


r 


it pt2 


\ 


-;) 


| 
Tee cA TT Uy + yA aide 


A ad (UE a TS S 1) 
ee otek jh, 


Spe 


etl “f 


| - x, | Uae (L-z+i)A | 
(3-72) n@reaa 


(pti-p) FP" A 


2 abes +ca] 


Sent 1) Cxe(pti-s)Aye 
eT (ptr-j) Piet! (p-k) | 


kr 


The only terms involving U, are 


- X, U3 
@ 


Ust (t-r)A 


and the inversion of this term yields: 


ay oe ea 
aa a eC : Oca 
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} 


The only term involving xg is in the previous inversion, neice. 


Le il £* dx, = 
rT 


we 


-~Tla~r\a P+! | = 
e ame) (2-peme > iene? 


Cea) remem!) 


@~ * (v-j4'\A 


pri 
f va 


Ce Pa!) 


-X2 (y=-pyA i p- 
oa CS [xi Paina) A] ; 
( Pei-gyr s Cz 0 CP=-\)! 


The first term involving Xo is 


em Xr (rm-yri\) 
and 
ae Lv -y+thA =7 (r-j+1\A 
S dx. — Ls A 
: (y-j +1NA 
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The next tenm is to note for K =p 


i -x, (r-e) A -T(r-p)A 
\ oe = » [i-e | 


lal (per-qyert ert ') Coaieg) PF! Q Pt? (pap) 


and for all other terms inside the sum 


Ik ep-k -Xq (r-p)A 
(=1} Ko eC 


P-! 
\ (oti) Kt! \ KT! | (e- kk) | 


=| 
I< 
(-1) a *"” 
Pa Cet lame ba lPlee (fap) ea 


- -p) A 
Lee . 
(ptt Ga) a a 
S eo 
i 
o (p-K- 1) 1 Clr-py a] | 


igacrefiore the probability thatyNa- = r iss 
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-T(t-rA ch cl df 5 at | 
ey | Cy 70 22) lek cone: 


J= p+2 
-tTlr-yai)a -T(r-p)d 
(4-pai) Ot NOt (yt 1) Cprr-s)PM APT (yay 


ae (-1)k ‘a =r p te EOL p-k-+ 
2 oar pag eta araeics kel ace pirat 


From the equations thus far derives the moments of the bivariate 


distribution can be developed as 


iB 
E(wrl= Dr PL ager] 
r=o0 
E{e,j=- 2&I gtser) | (48) 
OU eo 


a Fi Foe Ps) 


(=O 


4 -Us 
Ej By = ofa 





Us ~ 


NOW eel View Of tne convention tor UNGe e+ Caan), this 
> Oe ! 3 


bivariate density can be derived from equation C407, in that 


i =~ 
* 
trig (b= J dx Jo des Tus wus f (46 
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Where f* is the appropriate Laplace Transform depending on the 


value of m. To start let m< p and hence 


aE \ 
Sea | wet 
x 
f¥iuusus) = e cA » ee. 
Ure Us + on ast ca 


Where e USC! is the Laplace Transform of the degenerate 


random variable BT. ‘Therefore: 


m | 
When m = 0, 1, 2, ... p3 and b=cT 
Now for m = ptt 
~Ue (c-1)T Pri 
o (uy, 4s, Ug) = ‘ cA (c-N\A 
Wa + Use bor cn uy > (C-N)A 


= ~U (C= ae peu > mee yk Cust As + onl 
| nt i See SOT Cn ee oe en 
PI 


~uglc-1)T Pe! > =%e (ucts) — SNaie ix 
C e eC. 


Pr2 
(c-1) NN me 
Pt 


Tu,u,f* = G 


— 
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wl ay \ - Ugtdr = 2 
gle DT (. yy PAh Pe our Mae bs \ on Tc 


ifr dx, = ¢ 


T cO 
dx, ( dX, ale wy +* = ea eas ) (oN) . P 
ae 


(uy +A) P 
oe 2 cae > gan 
(Ug + )P*! foo WP)! (Used) f 
Inverting on us yields: 
Nae =a (p-nle b) = 
~ bd Pt ° 
27) te OT LS eG econ (48) 


P| 


It should be noted that another term would have appeared in this 
inversion except for limitation (c-1)T £ b & cT. 
By a development similar to that of equation (48), it can be 


shown: 
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tue. Be (™, Se 


ae S ener 


(A= myl Cy-p-2)! (m-y dV! 


-b> p 
e g 





Je Pre 





m -} w= 4 (49) 
Sey - [b-Tu en | 
Ci Rey Cay oe ae 
\ aa] \< 
bee T Kj -Aal > (-)) 
(m=) | er (p+ir-yyr' (m-pyP 8" 
Where m=pt+2, pt3, ..., i and all terms in Square brackets 


are followed by a unit step function with the same argument. 

The reliability program manager nowvhas all the distribution 
theory necessary in this case to assist him in making decisions 
about the testing program. The formulations look somewhat 
unmanageable, but with specific numbers these formulas should not 
be too difficult to solve. If the numbers become large, use of the 
normal approximation derived by Read [Ref.5] would probably be 


more practical. 
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APPENDIX A 


THEOREM ONE 


Theorem one is stated as follows: 


If f*(s) is the Laplace Transform of f(t) such that 


2 és 1 
= —- 
tS) (sta) stb) >a Fx b., them 


- bt eee V\-— | n-i-k 


f {+} = —_ , > en Lt (b-a) | (50) 
(a-b ) (b-a)™ (n-i1- k )1 


\<=0 


Proof: 


By the theorem of residues of complex variables, 


St 


e 
st fe 
CC ————— oo d\s+b ) (51) 
(sta)” cls (n-1) | ae a S=- Ol 


Now by the binomial expansion for derivatives of products, tf 


ast 
g(s)= —sp> 


then the (n-1)S* derivative with respect to s is: 





—| Sot n=! = - me 52 
g or an Di? (seb) FA a = 
nae 


Now substituting equation (52) into equation (51) yields: 
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st 
c 


= (sta) S=—-b + 
‘|. eo CO aa :. a ml Be 1)! 
ae aes j Jae 
Git Gry 2 en (n-1-j)!_[s=-a 
Hence: 
— Roe S ; 
= S. yo i=. 
¥ (azb)" (b-a) = j=0 r (i= ~ 
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Ae PeaNDIx B 
MAXIMUM LIKELIHOOD ESTIMATORS ie ONES i rare 


FAILURE RATE 


The estimator Ay for N upon completion of the test program 
G.€. Upon Obtainineipnaluliices) 71s developed as follows: 


The likelihood function, L( ) ), based on the samples of Sj> is: 


d ~CAS ~(C+p-yst)AS 
P-yri 
Lin) = J] oA ly Cespeyny Ae 2 
| yt yzpr! 
5 ul ys 
4) LOSS —~ (C+p-yt! f 
Be ees +e @ dtl iam (atp-yrt) dA © 


mire pri 


Now letting @(C } )=Ln¢« LC J )),» yields: 


> Y 
B(drA) = PbLn(ed) - coh¥ sy + > Lu { (eep-ys dA] 
j= Jape 
2 
— w (erpr-ysild Sy 
J pie 
P : 3 - 
/ ge | 
Bi (r\ = a -¢D § +> me -3 (C+P-ysl) Sy y 
Jt here j= Pt) 
Setting equation (53) equal to zero and solving for ) yields. 
° - 
an | 7 , 
ie € 7 33 i (ucaais y+ Sy 
rs Jer 
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and recognizing the right-hand side as B, yields the estimator for 


to be: 


XN 
A, = os (54) 


Br 


Now the maximum likelihood estimator for y\ at some fixed 
time T , prior to completion of the testing, where N 4 =Mép is 


given by: 


mM 
i 
By 


» > 


In the case wnere Nw =m > fp, the following data is 


available: T,, T Letting: 


vu 


>o- lo sane 


and noting that the probability of zero failures in(% Tali aaDss asa 
function of A > 1S: 


_ enim) A Ct fa) 


the likelihood function becomes: 


Pp aa d a bn 
mora y= KI civike 1 1 | Ghee 1) r. 
i j=ptt 


sie) A Si-E-.--)) ACT Ep 
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Now letting @ A ) = Ln(Lc¢ A )), yields: 


p mM ; 
Bor) = ple A )-oA 2, Sj } ee aitGtp=jbeA)) — 


an (tpt) G aS p- (CaGMERD helene 5) 
Je Obaat 


p mM mM 
J) - =e ane y S: + ? _l - > (Cis ine 
A =1 J alone iy, j=pt1 J 


Lo 
Lo. 


= (e-(m=-p)) (Y-T_) 


Setting @'( A )=0, yields: 


Mio 


aS 
A j 


Recognizing the right-hand side as B~ yields the estimator 


mm 
) Zot (ctp-jt1)S , tame -(IM—P)) CT —Tap 


for ~A to be: 


“~ 
nA 
a B37 f ’ (55) 
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